Ion-beam modification of the magnetic properties of GaMnAs epilayers 
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We study the controlled introduction of defects in GaMnAs by irradiating the samples with 
energetic ion beams, which modify the magnetic properties of the DMS. Our study focuses on 
the low-carrier-density regime, starting with as-grown GaMnAs films and decreasing even further 
the number of carriers, through a sequence of irradiation doses. We did a systematic study of 
magnetization as a function of temperature and of the irradiation ion dose. We also performed in- 
situ room temperature resistivity measurements as a function of the ion dose. We observe that both 
magnetic and transport properties of the samples can be experimentally manipulated by controlling 
the ion-beam parameters. For highly irradiated samples, the magnetic measurements indicate the 
formation of magnetic clusters together with a transition to an insulating state. The experimental 
data are compared with mean-field calculations for magnetization. The independent control of 
disorder and carrier density in the calculations allows further insight on the individual role of this 
two factors in the ion-beam-induced modification of GaMnAs. 

PACS numbers: 75.50.Pp, 75.60.-d, 74.62. Dh 
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I. INTRODUCTION 



During the last decade, III-V diluted magnetic semi- 
conductors (DMS) such as GaMnAs have been the sub- 
ject of intensive research, specially due to their pos- 
sible applications in semiconductor-based spintronics. 
Gai_a;Mn2,As presents long range magnetic order be- 
tween spaced Mn ions, mediated by holes [H-Q. The exis- 
tence of these holes in the band structure of the material 
is undoubtedly due to the incorporation of Mn ions. Con- 
sequently, crystalline defects have a crucial role in both 
electrical and magnetic properties of ferromagnetic semi- 
conductors and in device applications j^j . Continuous re- 
finement of growth techniques has led to a significant im- 
provement of the upper limit for the incorporation of Mn 
atoms at appropriate sites. This improvement towards 
the high-carrier-density regime has been very effective to 
enhance the Curie temperature of Gai_a;Mn2,As j5|]. On 
the other hand, some experiments studying and propos- 
ing prototypes of actual spintronic devices also point to 
the importance of the low-carrier-density regime [6|] . This 
is the case, for instance, of the proposed control of GaM- 
nAs magnetic anisotropy by the application of electric 
fields ^6,]. 

Different paths have been used to experimentally ac- 
cess the effects of carrier density and structural disorder 
while keeping the Mn concentration constant. Annealing- 
induced changes are very effective to enhance Curie tem- 
perature and are used to control magnetization and mag- 
netic anisotropy [H, 0, S]- However, it is not trivial to 
quantify the effect of a given annealing procedure on the 
density of defects. More recently, the addition of hydro- 
gen to the epilayer has been used to change GaMnAs 
magnetization properties @, In this latter case, it is 
possible to quantify the creation of defects by estimating 



the amount of hydrogen aggregated to GaMnAs. A third 
path is the implantation of heavy ions in the GaMnAs 
sample. Beams of Ga+ of 30 keV have been used to in- 
troduce deep trap levels in the epilayer [ll|. In the two 
latter methods additional material is incorporated in the 
GaMnAs sample. We use a different approach to study 
structural defects in GaMnAs films. 

We produce defects in the GaMnAs sample in a con- 
trolled way by irradiating the samples with ion beams 
fast enough to cross the GaMnAs layer. Thus, the pro- 
jectile ions are implanted only in the non-magnetic GaAs 
substract. However, going through the GaMnAs film, 
the fast ion beam leaves behind defects in the crystalline 
structure, enhancing the disorder of the system. The 
same technique has been applied in the study of radiation 
tolerance of semiconductors and devices that are used in 
solar cells, detectors and satellites. Consequently, there is 
a great amount of information on the characterization of 
defects produced by energetic ion beams. They have also 
been used in different contexts: the study of the effects 
of structural disorder in superconductors [HI-Gll, ferro- 
magnets [l6l |. n-type and p-type semiconductors. 

In this article, we study the manipulation of carrier 
density and structural disorder of GaMnAs by the use 
of MeV ion beams. We focus on the low-carrier-density 
regime, decreasing this density as defects are introduced 
by the ion beam. A second effect of the ion-beam irradi- 
ation is the increase of structural disorder controlled by 
the irradiation dose. We present magnetic and transport 
measurements of GaMnAs epilayers irradiated with ener- 
getic ion beams and analyze the effect of the irradiation 
on the magnetic and electric properties of the DMS. 

The article is organized as follows: in section [IT] we dis- 
cuss the irradiation technique, the defects normally cre- 
ated in semiconductors by this process and the method 



2 



that can be used to quantify the formation of defects in 
the DMS due to the irradiation. X-ray diffraction mea- 
surements are used to discuss the enhancement of the 
density of defects in both the epilayer and in the sub- 
strate. In section mil we present room temperature sheet 
resistance measurements as a function of the ion dose 
and discuss the similarities between GaMnAs and other 
p-doped GaAs semiconductors. In section ITVl we present 
the magnetic measurements performed in samples with 
different irradiation doses and discuss how the irradiation 
process affects the magnetic order in GaMnAs. Finally, 
in section |V] we use a mean-field approach for an impu- 
rity band model to compare our experimental results with 
the theoretical magnetization curves and understand the 
roles of the decrease of carrier concentration and of the 
enhancement of the structural disorder in the modifica- 
tion of the magnetic properties of the DMS. 



II. IRRADIATION PROCESS 

In n-type and p-type GaAs semiconductors, the beams 
introduce a rich variety of defects. The defects in n- 
type GaAs are better understood than the ones created 
in p-type GaAs. However, it is well established that in 
both cases these defects in GaAs semiconductors reside 
in the As sublattice and most of them are primary de- 
fects (related to vacancies and interstitials) ^Igj. The 
irradiation process produces similar quantities of defects 
in the Ga and As sublattices but the vacancies in the 
Ga sublattice tend to recombine immediately with the 
interstitials, since they have opposite charge. We do not 
have any previous information on the role of Manganese 
atoms in the formation of complex defects in irradiated 
GaMnAs. However, previous studies of such defects in 
non-irradiated samples suggest the possibility of forma- 
tion of isolated Mn interstitials with As nearest neigh- 
bors and pairs of Mn interstitials with As nearest neigh- 
bors [la,^]. 

The range of the implanted-ion profile is determined 
by the choices of projectile atomic number (Z) and ki- 
netic energy. This range is independent of the beam dose 
(measured in ions/cm^). On the other hand, in order to 
determine the density of the ion-induced defects in the 
relevant sample epilayer, knowledge of the dose is also 
necessary. Ion current density and irradiation time are 
the parameters used to experimentally control the beam 
dose. 

The Gai_a;Mna;As samples are grown on a GaAs sub- 
strate by Low- Temperature Molecular Beam Epitaxy 
(LT-MBE). The epilayers have 200nm or 100 nm thick- 
ness and a Mn concentration, x, of 5%. Hall mea- 
surements performed with non-irradiated samples deter- 
mined a density of holes of 7x10^^ cm^"^ for the 200nm 
and 2x10^^ cm~'^ for the 100 nm samples. We irradi- 
ate the samples, at room temperature, with ion beams 
from a 1.7 MeV NEC Tandem Accelerator. As the Z of 
the used ion increases, the dose necessary to create the 



same density of defects decreases. For decreasing ener- 
gies, there is an increase in the creation of defects. In 
order to cover a large range of defect densities, we use 
100 keV protons (4 different doses), 1000 keV protons (3 
different doses) and 700 keV Li+ ions (4 different doses). 
Ion-beam current densities range from 0.21 pA/cm^ to 
1 nA/cm^ and ion-beam doses are in the range between 
1.3x10^ ions/cm-^ and 6.7x10^'' ions/cm^. 

Both the ion beam implantation profile and the va- 
cancy density profile are simulated using the SRIM 2008 
code [2l| . The complexity of the cascade collisions in the 
solid target limits the accuracy of the simulation for de- 
fect production. However, reasonable agreement between 
SRIM simulation and experiment has been obtained for 
the vacancy production in GaAs [22j. We used 15 eV for 
the displacement parameter in the SRIM code. Although 
there are no previous simulations for GaMnAs, this value 
is well accepted for GaAs simulations [23] . 

Figure [T] shows the SRIM simulation results for 700 
keV Lithium projectiles. It shows that almost none of 
the projectiles are implanted into the GaMnAs. The in- 
set at Fig. [T] indicates a relatively uniform defect creation 
through the epilayer. The integral of each curve of the in- 
set in Fig.[T]gives the number of vacancies per ion created 
in the epilayer for the correspondent target component. 
These numbers, multiplied by the beam dose and divided 
by the epilayer thickness, give the average volume den- 
sities of produced crystal vacancies for Ga, Mn, and As. 
We use these densities to quantify the induced disorder 
in the GaMnAs epilayer. Although SRIM has to rely 
on theoretical simulations of the motion of atoms in the 
irradiated sample, this procedure still gives a quantita- 
tive scale for the introduced disorder that temperature 
changes can not provide. 

It is important to point out that SRIM simulations do 
not consider temperature and do not account for recombi- 
nation processes. This can be clearly seen in Fig.[l] where 
the simulation shows a considerable number of Ga vacan- 
cies created by the irradiation, although we know that 
most of these defects recombine. More specific molec- 
ular dynamics simulations for irradiated GaMnAs, like 
those reported for GaAlAs [l^, would be desirable but 
are not available in the literature. However, the effect of 
this recombination is roughly dose independent, except 
for extremely large doses that can cause amorphization 
even in the GaMnAs layer. Thus, recombination pro- 
cesses would basically change our chosen defect scale by 
a constant factor. 



A. X-ray diffraction of ion-beam irradiated samples 

X-Ray Diffraction (XRD) is routinely used to charac- 
terize MBE-grown Gai_:cMna;As epilayers (e.g. (23 - [26t 
). These measurements provide estimates for changes 
in lattice parameters relative to bulk GaAs and, there- 
fore, for Mn content. Furthermore, XRD is also a tool 
that has been used to experimentally access the effects 
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FIG. 1: The depth profile of 700 keV Li+ ions implanted in 
the target. A negligible fraction of incident ions is implanted 
in the magnetic GaMnAs epilayer. The inset shows the num- 
ber of created vacancies / Li"*" ion / nm, for each target ele- 
ment (Ga, Mn, and As), as a function of the epilayer depth. 
Simulations obtained with the SRIM 2008 code. 



of ion-beam irradiation on bulk GaAs samples [27|, |28j . 
In this latter case, XRD measurements and data analy- 
sis based on dynamical diffraction theory result in infor- 
mation on the related depth profiles for implanted ions, 
induced structural defects, and strain. Our irradiated 
samples combine characteristics of both kinds of sam- 
ples mentioned above: non-irradiated GaMnAs samples 
and bulk implanted GaAs. Thus, it is convenient to use 
XRD measurements to pave the way for the transport 
and magnetization studies of irradiated GaMnAs epilay- 
ers presented in the following sections of this paper. 

Figure [5] shows X Ray Diffraction Rocking Curves for 
pristine and for irradiated 100 nm GaMnAs epilayers 
grown on top of GaAs. The data was taken by a Bede 
200 x-ray diffractometer equipped with an Enraf Nonius 
x-ray generator. A copper tube source was selected, and 
the x-ray Ka line was filtered using a Ni foil. Like in all our 
samples, there is also a 100 nm LT-GaAs layer between 
the GaMnAs layer and the GaAs Bulk. The irradiation 
dose is 1.5x10^^ Li ions/cm^ and the Lithium beam en- 
ergy is 700 keV. The pristine sample shows in Fig. [2] the 
typical behavior of MBE-grown Gai-ajMn^jAs epilayers. 
There are two main structures. The larger and narrower 
peak corresponds to the GaAs substrate (with a merged 
contribution from the 10 nm LT-GaAs layer [23|) and a 
smaller and wider peak corresponding to the GaMnAs 
epilayer. The other structures are interference fringes. 
Those fringes can actually be used to estimate the width 
of the GaMnAs epilayer (100 nm in this case). 

Figure [2] also shows an additional set of peaks charac- 
teristic of Rocking Curves for ion-beam irradiated GaAs 
bulk samples. Deep into the sample, far from the GaM- 
nAs epilayer, Li ions are implanted (see Fig. [1} leading 
to an approximately gaussian depth profile of cristaline 
defects. This results in a continuous modification of the 
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FIG. 2: Comparison of measured 2Q/Q scans for symmetric 
(004) reflex of 100 nm Gai_2,Mna; As epilayers, with x—0.05, 
grown on GaAs. The angle 2Q = 66.048° corresponds to 
the (004) reflex of the GaAs substrate. The curve with a 
narrow GaAs bulk peak corresponds to pristine sample. The 
curve with a set of peaks for GaAs correspond to a sample 
irradiated with 700 keV Li^ ion beam, with a dose of 1.5x 10^^ 
Li ions/cm^. The inset shows a comparison between our data 
for 700 keV Li^ and data for a projectile with equivalent 
implantation range, 2 MeV ions, incidend on GaAs bulk 
sample [2^ . 



GaAs lattice parameter |27l . |28| with a maximum strain 
modification around 2000 nm. The boundaries of this 
buried layer of defects are not as sharp as in the case of 
the MBE-grown Gai_a;Mna;As. However, the implanted 
layer of defects is well enough defined to result in a sec- 
ond set of fringes of interference (see Fig. [5]). The ef- 
fective width of this layer can also be estimated from 
the distance between the fringes and is approximately 
650 nm. This value is consistent with the implantation 
depth profile shown in Fig. [T^nd obtained from SRIM 
simulation. In the inset of Fig. 2 we compare, in a linear 
scale, Rocking Curves for our irradiated sample and for 
bulk GaAs irradiated by 2 MeV oxygen ions The 
700 keV Li"*" ions and the 2 MeV 0+ ions have approxi- 
mately the same implantation range and result in similar 
bulk Rocking Curves. It is important to note, however, 
that the amount of defects (and therefore strain) intro- 
duced near the surface is much smaller than in the bulk. 
The previous study of Xiong et al with O'^ beam on bulk 
GaAs ^28] also shows that, even for depths around the 
maximum of deposition the Rocking curves shown at the 
inset of Fig. [2] correspond to only a moderate dose in 
terms of the density of defects created. An implantation 
dose 200 times larger would be neces sary to create an 
amorphized layer of GaAs in the bulk [28| . 

Regarding the GaMnAs peak, no change in position or 
witdth, within experimental errors, is noticed from pris- 
tine to irradiated sample in Fig. [2j The unaltered diffrac- 
tion angle shows negligible changes in lattice parameters 
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of the GaMnAs layer. Previous studies of ion-beam irra- 
diated 2 mm GaN epilayers [29|] show a smaU increase of 
the peak width attributed to the enhancement of defect 
density. In our irradiated sample such a modification 
is not observed. In summary, the XRD measurements 
show that the ion-beam irradiation docs not change ap- 
preciably the lattice constants of our GaMnAs epilayers. 
The passage of the ion beam through the superficial layer 
moves part of the atoms from their initial positions. A 
large variety of defects is created. The average rearrange- 
ment, however, does not result in a relevant change of the 
Gai_a;Mn2,As crystalline structure. On the other hand, 
this rearrangement does modify transport and magnetic 
properties of the epilayer. Particular defects, like the 
transformation of a substitutional Mn in an interstitial 
Mn, are certainly more decisive than others in the mod- 
ification of these properties. In the next two sections 
we discuss measurements that can access those ion-beam 
induced changes, and the correspondent experimental- 
technique limitations in this specific case. 
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FIG. 3: Sheet resistance (SR) measured as a function of the 
700 keV Li^*" ion dose for a 100 nm thick GaMnAs sample. 



III. IN-SITU SHEET RESISTANCE 
MEASUREMENTS 

To facilitate the comparison between the irradiation 
process on GaMnAs and results on p-GaAs, we per- 
formed sheet resistance Ks measurements. The i?^ mea- 
surements are routinely used to characterize doped non- 
magnetic semiconductors modified by li ght energetic ion 
beams in the so-called isolation process jSOl [Slj . In that 
case, the use of ion beams induce a dramatic decrease of 
carrier concentration in selected area sectors of the epi- 
layer. This procedure keeps the planarity of the sample 
and is used as a tool in the microelectronics industry. 

In order to calibrate our equipment, we performed a 
first irradiation of a GaAs sample (5-doped with C. We 
compared the doses necessary to electrically isolate the 
sample with the ones of Ref. [s^ and their results were 
reproduced. Subsequently we performed the irradiation 
of the DMS samples with different ions and energies. 
Combining the different projectiles we span 6 orders of 
magnitude in the introduced defect density. The irra- 
diation by 1000 keV protons with doses up to 6.7x10^^ 
ions/cm^ did not produce any considerable change in the 
properties of the DMS. High doses of 100 keV proton 
produced a very small change. The choice of Li"*", an ion 
with larger atomic number and mass, finally produced 
visible changes in the magnetic and transport properties 
of the samples. 

In this experiment, the Rg values were measured in- 
situ after every irradiation step of the dose accumulation. 
The measurements were performed at room temperature 
using the van der Pauw method |33|. Figure [3] shows the 
sheet resistance as a function of the 700 keV Li+ ion dose 
for a 100 nm thick GaMnAs sample. Different behaviors 
can be identified according to the range of doses analyzed. 

At very low doses Rs increases slowly. As the dose in- 



creases, there is a pronounced enhancement of the R^ up 
to a maximum and then a decrease with a power law of 
the dose. This behavior is also found in non -magnetic 
epilayers of both n and pdoped semiconductors grown on 
top of undoped wafers [3l|. Bellow the maximum of Fig. 
[21 the increase of the Rg is a direct consequence of the 
decrease of the density of carriers in the epilayer caused 
by crystalline defects introduced by the trespassing ion 
beam. To understand the high-dose regime, however, it 
is necessary to take into account what happens to the 
bulk. The Li ions are buried in the GaAs substrate at a 
depth of approximately 2000 nm (see Fig. [Ij . For high 
ion doses they create a deep layer of modified GaAs in 
the bulk with a concentration of defects so high that the 
conduction by hopping between GaAs defects becomes 
relevant. Therefore, a simple model (e.g. Ref. [11]) con- 
sidering two resistances in parallel (the GaMnAs sheet at 
surface and the high-defect-density GaAs sheet deep at 
2000 nm) can explain the dose dependence seen in Fig. [31 
It is important to note one peculiarity of the effect 
of incident ion beams in the case of the GaMnAs epi- 
layer: the buried implanted layer has no magnetic prop- 
erties. Thus, magnetic measurements, which infer the 
hole-mediated ferromagnetism are expected to refiect a 
smooth decrease in carrier densities coming solely from 
the GaMnAs epilayer. On the other hand, additional 
effects, as the increase of structural disorder, may be in- 
duced by the ion beams and compete with carrier con- 
centration changes to modify the magnetic properties of 
the DMS. 



IV. MAGNETIZATION MEASUREMENTS 

For the magnetization measurements discussed in this 
section, we used a 200 nm thick sample with a higher 
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concentration of carriers than the one in section IIIII 
but the same nominal concentration of Mn. After each 
ion irradiation, we perform magnetic measurements in a 
SQUID magnetometer. In order to analyze the sample 
anisotropy, an in-plane magnetic field is applied on the 
sample with two orthogonal orientations: parallel and 
perpendicular to easy axis ([110]) of the non- irradiated 
sample. 

In Fig. m we show the value of the [110] magnetiza- 
tion at 10 K as a function of the number of created Mn 
vacancies obtained by the SRIM simulation, which also 
measures the amount of disorder in the system. Here, 
we clearly see a universal curve, which is independent of 
the type or energy of the ion beam. The total density of 
Mn vacancies is used here as a scale. The SRIM simu- 
lations do not take into account recombination processes 
and these numbers can be considerably smaller than the 
ones shown here. 

Figure [5] shows the easy-axis magnetization versus tem- 
perature measurements for the samples irradiated with 
Li+. We can clearly see a decrease in the magnetiza- 
tion at low temperatures, even though the critical tem- 
perature is only slightly modified. Up to now, most of 
the discussion was focused on the consequences of the 
variation of resistivity due to the irradiation process but 
we can also understand the role of the increase of disor- 
der in GaMnAs. For an increasing amount of disorder, 
we observe a change in the concavity of the magnetiza- 
tion curves (see FiglS]), which is compatible with highly 
compensated GaMnAs samples. The suppression of the 
magnetization for low T is not followed by considerable 
variation of Tc- Previous theoretical calculations predict 
that for a fixed number of carriers, an increase of disor- 
der will lead to a decrease of the magnetization at low 
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FIG. 5: Magnetization versus temperature for samples irradi- 
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of the non-irradiated sample at 10 K. The figure shows that 
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TABLE I: Properties of GaMnAs sample irradiated with 700 
keV Li^ ions: created Mn "Vacancies (SRIM simulation) and 
sheet resistance (at room temperature) 
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temperatures and a considerable increase of Tc |34l . |35j . 
Here, any variation of the amount of disorder is accom- 
panied by a change of the resistivity, as can be seen in 
Table U and consequently, the carrier concentration. As a 
decrease of the carrier concentration decreases Tc, these 
two effects compensate each other and Tc remains almost 
unaltered until very high doses of ion beams. 

In Fig. [21 for the most intense irradiation (D) we see 
an unusual decrease of the magnetization for low tem- 
peratures. In principle, this feature could be attributed 
to the creation of shallow levels in the DMS where the 
trapped carriers can be thermally activated 36] . Another 
possibility is the existence of a blocking temperature due 
to the formation of magnetic clusters in highly irradi- 
ated samples. To further investigate these possibilities, 
we performed standard zero-field cooled (ZFC) and field- 
cooled (FC) magnetic measurements. Fig. [5] shows that 
the ZFC and FC measurements have different behaviors 
at low temperature. These results rule out the picture 
of an increase in the density of holes due to thermally 
activated carriers. The existence of a blocking temper- 
ature Tb that depends on the applied field is consistent 
with the existence of magnetic clusters. After the irradi- 
ation process, due to the reduced number of carriers and 
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disorder, the holes are locahzed, similarly with what is 
observed for GaMnAs films at low Mn concentration [sl] . 
Consequently, the DMS is composed of weakly interact- 
ing clusters, probably close to the percolation transition. 
For low temperatures, the ZFC measurements can freeze 
the system in a random orientation of these ordered mag- 
netic clusters. The system needs to overcome a barrier 
of ksTs ~ 0.1 meV in order to reorientate the mag- 
netic clusters, as illustrated in the inset of Fig. ^ Simi- 
lar behaviour has been predicted by Monte-Carlo simu- 
lations [st! ]. 

We can further analyze the effect of the irradiation 
process by estimating the fraction of Mn atoms that are 
participating in the magnetism. The effective value Xef / 
in Gai_a;Mna;As is the concentration of Mn atoms con- 
tributing to the magnetic order. If we make the assump- 
tion that the interstitial Mn are paramagnetic and do not 
contribute to the measured magnetic moment, we can es- 
timate Xef f from the saturation magnetization using 
Ms = ixeffgfJ-BS/Vf^QiY where g is the gyromagnetic ra- 
tio, /iB is the Bohr magneton, V^^qYI is the volume of the 
unit cell and 5 = 5/2 is the spin of Mn. It is impor- 
tant to notice that here we consider the concentration of 
of positional Mn as the Xef / • For high nominal concen- 
trations of Mn, this assumption is less accurate, as part 
of the interstitial Mn are coupled antiferromagnetically 
and Xeff is smaller than the concentration of positional 
Mn [ig. 

Figure[71[a) shows values of Xeff obtained from magne- 
tization measurements as a function of the dose of inci- 
dent 700 Kev Li"*" ions. The irradiation, on average, pro- 
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FIG. 7: Parameters experimentally determined from magne- 
tization measurements, as a function of the incident 700 keV 
Li"*" ion dose: (a) the effective value of x in Gai_a;Mn2,As (b) 
the density of Mn vacancies created by the ion beam. Due 
to the formation of magnetic clusters in the highly-irradiated 
sample (see Fig. (6]), the estimate of the effective concentration 
for this particular experimental point is less accurate and the 
quantitative estimate is less reUable. 
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FIG. 6: (a) ZFC (closed symbols) and FG (open symbols) 
magnetization x temperature measurements for two different 
applied fields: 50 Oe and 200 Oe. Inset: Schematic repre- 
sentation of the blocking temperature barrier separating a 
configuration of random oriented magnetic clusters from a 
configuration of magnetically oriented clusters 



gressively removes Mn atoms from sites where they are 
active for ferromagnetism. Although the nominal value 
of X for this sample is 5 %, the Xeff before irradiation is 
only 1.2 %, a difference that can expected for as-grown 
samples [4^. For the highest dose irradiation, Xeff de- 
creases by about one order of magnitude. However, the 
estimation of the effective concentration in this sample 
is less accurate once Fig. [5] shows different values of the 
magnetization at low temperatures for ZFC and FC mea- 
surements. It is possible that for the saturation magneti- 
zation we still have isolated magnetic clusters with differ- 
ent magnetic moment orientations. Although more Mn 
participate in the magnetism, this is not reflected in the 
total magnetization of the system. Magnetic frustration 
is also a source of the underestimation of x^ff 41]. 

In order to compare these experimental results with 
computational simulations like those using the SRIM 
code [2ll |. it is convenient to present the data in an al- 
ternative way. Figure [7] (b) presents the density of active 
Mn atoms removed from their initial sites as a function 
of the dose of incident ions. The derivative of this curve 
provides the rate of introduction of defects, Ti, (i.e., num- 
ber of created Mn vacancies per incident ion per unit of 
length). This quantity is also an output from SRIM cal- 
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culations and is sometimes used as an input parameter for 
theoretical modeling of the effects of ion beams in elec- 
trical properties of solids. This modeling often assumes 
that the variation of concentration of a certain kind of 
defect depends linearly on the implantation dose. 

Our defect-specific experimental data (for Mn vacancy 
creation) offer a way to benchmark calculations of Ti and 
to test the above-mentioned linearity assumption. In Fig. 
[TJb) the linear regime can be seen for doses up to ap- 
proximately 1.5 xlO^'' Li+ ions/cm^. This behavior ob- 
viously has to change for higher doses because otherwise 
the number of displaced Mn atoms would increase be- 
yond the total number of Mn atoms. Our measurement 
at 6.7 xlO^'* Li+ ions/cm^ is clearly beyond the linear 
regime. Close to the origin, corresponding to Xe//=1.2 
%, we experimentally obtain Ti= 6.1 x 10~^ Mn vac/(ion 
angstrom). The SRIM calculation, assuming displace- 
ment energy of 15 eV for all atoms, give us Ti= 8 x 10~^ 
Mn vac/ (ion angstrom). The SRIM code simulations and 
our estimates therefore agree within 30%. 



V. HOLE CONCENTRATION x STRUCTURAL 
DISORDER 

An effective approach to understand the effect of the 
irradiation process is to compare our experimental results 
with theoretical calculations that allow for the variation 
of carriers and spin concentration and also take into ac- 
count the structural disorder. Since we are dealing with 
small carrier and Mn concentrations in as-grown sam- 
ples, an imp urity band model is better suited for this 
analysis [42l444l | . The Hamiltonian is given by 



situations of disorder: for the non-disordered case, the 
Mn ions form an ordered lattice. In the weak disorder 
case we allow the Mn ions to move to a nearest neighbor 
with a probability of 20%. Finally, in the strong disor- 
der limit, they are randomly distributed in the DMS. For 
comparison with the experiments, we use two concentra- 
tions, x=0.009 and x=0.004, where due to the compensa- 
tion mechanisms, the hole density Uh is a fraction of the 
effective Manganese density Hmu- For all the mean-field 
calculations we use tridimensional lattices with 1000 Mn 
spins and for the disordered systems we perform avarages 
over 500 realizations of disorder. 

We consider two different fractions, nh/riMn = 0.2 and 
nh/nMn = 0.25. From the experimental results, we ob- 
serve that although the saturation magnetization and the 
shape of the magnetization curve are strongly affected by 
the irradiation process, the critical temperature is only 
slightly modified (see for example, samples A and C). 
Using this mean-field theory we can see, as illustrated in 
Fig. [8] (a) , that if we decrease the Mn concentration with 
a fixed ratio nh/riMn = 0.25, we obtain both a suppres- 
sion in Tc and a reduction in the saturation magnetiza- 
tion without a change in the shape of the magnetization 
curve. On the other hand, if we keep x fixed and change 
the amount of disorder in the system, we can increase 
Tc, as is clearly seen in Fig. H] (b) for x = 0.009. The 



H 



where c]^ is the creation operator of a hole with spin a 
in the bound state associated with the zth Mn impurity. 
Ri {i = l,Nd) are the positions of the Mn impurities, 
and the hopping matrix = t{\Ri — Rj\) is given by 
t(r) = 2(1-1- r/as) exp {—r/aB)Eb fi^, where the Et is 
the binding energy of the hole, and as = efi^ /nihe'^ is the 
hydrogenic Bohr radius. The second term is the antifer- 
romagnetic exchange interaction between Mn spins S{i) 
and hole spins where Si = clg^\SapCip . The exchange 

integral is given by Jij — J exp (— 2|i?j — Rjl/as), which 
is related to the probability of finding the hole in the im- 
purity state around j on the zth Mn spin. The last two 
terms in Eq. ([T|) describe interactions of the spins of the 
Mn and holes with an external magnetic field H. 

We perform a mean-field calculation, following the 
same procedure described in references [33, Eal with car- 
rier and Mn concentrations similar to the experimental 
data. In our calculations, we use as — 7.8A, m/i — 
0.5me, Eb ~ 112 meV, J = 15 meV and the lattice con- 
stant of GaAs a = 5.65A. We consider three different 
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FIG. 8: (a) Mean-Field Magnetization versus temperature 
for X = 0.009 (solid line) and x = 0.004 (dotted line) in the 
absence of disorder, (b) Magnetization versus temperature 
for X = 0.009 and non-disordered (solid line), weakly dis- 
ordered (up triangles) and strongly disordered (circles) sys- 
tems, (c) Magnetization versus temperature for x = 0.004 
non-disordered (dotted line), weakly disordered (stars) and 
strongly disordered (down triangles) systems, (d) Magne- 
tization versus temperature for non-disordered x = 0.009 
(solid line), weakly disordered x — 0.009 (up triangles) and 
strongly disordered x — 0.004 (down triangles). All curves 
use Uh/nMn = 0.25 and are normalized by the magnetization 
of X = 0.009 at zero temperature. 
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FIG. 9: (a) Mean-Field Magnetization versus temperature in 
the absence of disorder for: (a) x — 0.009 with Uh/riMn = 0.25 
(soUd Une), nh/riMn = 0.15 (dashed hne) and nh/riMn = 0.10 
(dotted hne) and (b) x — 0.004 with rih/riMn = 0.25 (sohd 
hne), rih/riMn = 0.15 (dashed hne) and nh/riMn = 0.10. AU 
curves are normahzed by the magnetization of a; = 0.009 at 
zero temperature. 



disorder can also modify the shape of the magnetization 
curve, as we see in Fig. ^c) for x — 0.004. So if instead 
of just changing the Mn concentration we also increase 
the disorder for the lower concentration, the change in x 
is compensated by an increase in Tc due to the increase 
in the disorder . We can also see a change in the con- 
cavity of the curve, in agreement with the experimental 
data. In this sense, in Figure [SJd) we emulate the irra- 
diation process by a suppression in the Manganese con- 
centration and an increase in the disorder of the system. 
The three magnetization curves are similar to the ones 
obtained experimentally for non-irradiated samples, low 
and high ion doses and can be compared with the results 
shown in Fig. [5l If instead of increasing the disorder of 
the system, we consider an increase in rih / umu after the 
irradiation process, we can also compensate the change in 
Tc- However, we do not obtain the correct magnetization 
curve shape. We could also consider that the irradiation 
process does not modify the Mn concentration but only 
decrease the carrier concentration but in this case we can 
understand the decrease in the critical temperature but 
not the strong suppression of the saturation magnetiza- 
tion, which does not depend on the carrier concentration 



at low temperatures. In Fig. [5] we support our analysis 
by showing the magnetization curves for both x = 0.009 
and X = 0.004 and different carrier concentrations. 
VI. CONCLUSIONS 

We irradiated GaMnAs epilayers by 100-1000 KeV 
light ion beams and studied the change in their trans- 
port and magnetic properties as a function of the ir- 
radiation dose. We performed in-situ sheet resistance 
measurements and found that irradiation by ion beams 
can electrically isolate a diluted magnetic semiconductor, 
similarly to what is seen in p and n-doped GaAs. We per- 
formed magnetic measurements on the irradiated samples 
and observed a suppression of the saturation magnetiza- 
tion with the increase of the irradiation dose. However, 
this same change was not observed in the critical tem- 
perature, which was only slightly modified. In order to 
better understand the roles of the change in both hole 
concentration and disorder of the system in the observed 
magnetization curves, we compared our experimental re- 
sults with a mean-field calculation for the impurity band 
picture. 

This comparison suggests that the increase of disor- 
der plays an important role in the magnetic properties 
of the DMS after irradiation. Furthermore, our experi- 
mental results show that the increase of disorder due to 
an intense irradiation process leads also to the forma- 
tion of magnetic clusters in the sample. We performed 
field-cooled and zero-field cooled magnetic measurements 
and found the energy barrier necessary to overcome such 
configuration. 

We show that the irradiation of GaMnAs epilayers by 
100-1000 KeV light ion beams can modify the sample in 
a quantitatively controlled way. This proof-of-principle 
experiment has interesting potential consequences for the 
manipulation of GaMnAs. It is important to note that, 
with an appropriate experimental setup, keV and MeV 
ion beams can be focused to few-microns-wide spots. 
This opens a door to, keeping the planarity of the sam- 
ple, writing on the GaMnAs and creating neighboring 
regions with different carrier densities. It is actually pos- 
sible, controlling the ion beam intensity and irradiation 
time during the writing process, to create regions on the 
sample where the carrier density varies continuously. 
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